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Abstract 

Critical phenomena and Goldstone mode effects in spin models with 0(n) rota- 
tional symmetry are considered. Starting with the Goldstone mode singularities in 
the XY and 0(4) models, we briefly review different theoretical concepts as well as 
state-of-the art Monte Carlo simulation results. They support recent results of the 
GFD (grouping of Feynman diagrams) theory, stating that these singularities are de- 
scribed by certain nontrivial exponents, which differ from those predicted earlier by 
perturbative treatments. Furthermore, we present the recent Monte Carlo simulation 
results of the three-dimensional Ising model for very large lattices with linear sizes up 
to L = 1536. These results are obtained, using a parallel OpenMP implementation 
of the Wolff single cluster algorithm. The finite-size scaling analysis of the critical 
exponent rj, assuming the usually accepted correction-to-scaling exponent u w 0.8, 
shows that 77 is likely to be somewhat larger than the value 0.0335 ± 0.0025 of the 
perturbative renormalization group (RG) theory. Moreover, we have found that the 
actual data can be well described by different critical exponents: r\ = u> = 1/8 and 
v = 2/3, found within the GFD theory. 

PACS 05.10.Ln; 75.10.Hk; 05.50.+q 

1 Introduction 

Critical phenomena are observed in vicinity of phase transition points in a variety of sys- 
tems (e. g., solids - ferromagnets, ferroelectrics; fluids - superfluid A-transition; liquid-gas 
transition, etc.), which manifest themselves in power-like singularities of physical observ- 
ables, described by critical exponents. In many systems, like isotropic ferromagnets (e. g., 
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polycrystalline Ni), certain singularity is observed not only at the critical temperature 
T = T c , but also at vanishing external field h — >■ for any temperature T < T c . This 
phenomenon, known as the Goldstone mode singularity, is also described by power-like 
divergences and certain exponents. In this sense, it can be viewed as a critical phenomenon, 
which takes place in vicinity of the critical line h = in the T-h plane. 

The vicinity of critical point (or line) is not the natural domain of validity of any 
perturbation theory, therefore one should resort the non-perturbative methods such as 

• Exact and rigorous analytical solution methods (transfer matrix methods, combina- 
torial methods, Bethe-ansatz), 

• Conformal field theory analysis, 

• Non-perturbative renormalization group (RG) analysis, 

• Numerical transfer-matrix calculations, 

• Molecular dynamics simulations, 

• Monte Carlo (MC) simulations. 

The exact solutions [HEIE] and conformal field analysis [I] are powerful tools to determine 
the critical exponents of two-dimensional models. However, these approaches are not 
helpful in most of the three-dimensional cases. The MC method can be used here. The 
exact non-perturbative RG equations are known for various models on fractal lattices (see, 
e. g. [5] and references therein), and also in four dimensions a rigorous RG analysis has been 
made (see e. g., [6]). In most of other cases, the non-perturbative RG equations provide 
only approximations, the accuracy of which cannot be well controlled (see, e. g., [7J El [9]). 
To the contrary, the MC method allows, in principle, to reach any desired accuracy. 
Besides, it allows to treat much larger systems than the other numerical methods. It 
is extremely important in studying critical phenomena. Therefore, the MC method has 
no real alternative in non-perturbative determination of critical exponents of the three- 
dimensional systems like lattice spin (0(n)) models. 

The general framework of perturbative methods includes the study of critical point 
phenomena at T — > T c , as well as the Goldstone mode effects at T < T c . Our aim is 
to verify nonperturbatively the validity of different theoretical perturbative approaches. 
Hence, the MC test of the Goldstone mode singularities is one of the important points 
here. We start with the Goldstone mode effects, since this case is simpler in the sense that 
it does not require the determination of the critical temperature T c , i. e., there are less 
ambiguities as compared to the T — > T c case. 

2 Goldstone mode singularities in the 0(n) models 

Here we consider such lattice spin models, where the spin is an n-component vector with 
n > 2. These are called 0(n) models due to the 0(n) global rotational symmetry exhibited 
by the n-vector model in absence of the external field. The Hamiltonian T~L of such a model 
reads 




(1) 
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where T is temperature, measured in energy units, Sj is the spin variable (n-component 
vector of unit length) of the ith lattice site, f3 is the coupling constant, and h is the 
external field with magnitude | h |= h. The summation takes place over all pairs (ij) of 
the nearest neighbors in the lattice. 

In the thermodynamic limit L — > oo below the critical point, i. e., at (3 > (3 C , the 
magnetization M(h) and the Fourier-transformed transverse (Gj_(k)) and longitudinal 
(G||(k)) two-point correlation functions exhibit Goldstone mode power-law singularities: 

M(h) - M(+0) oc h p at h ->> , (2) 
G±(k) oc AT A± at h = +0 and k , (3) 
G||(k) oc fc -A ll at h = +0 and k -> . (4) 

According to the standard theory [TQl [III H21 H31 El , A j. = 2 and Ay = 4 - d hold for 
2 < d < 4, and p = 1/2 is true in three dimensions. More nontrivial universal values are 
expected according to [15], such that 

d/2 < X± < 2 , (5) 
Ay = 2A± - d , (6) 
p = (d/X ± ) - 1 (7) 



hold for 2 < d < 4. These relations have been obtained in |15j by analyzing self-consistent 
diagram equations for the correlation functions without cutting the perturbation series. 
This approach is based on certain grouping of Feynman diagrams introduced in [16J, and 
therefore is called the GFD theory. 

MC simulations have been performed |17 } 118 } [19] to test these relations for the 0(2) 
and 0(4) models on the simple cubic lattice. The theoretical predictions ([5]) and ([6]) 
are confirmed by MC simulation results in the 3D 0(4) model [19J, where an estimate 
Aj_ = 1.955±0.020 has been found for the transverse correlation function. Moreover, it has 
been stated that the behavior of the longitudinal correlation function is well consistent with 
Ay being about 0.9, in agreement with (|6|) at A^ being about 1.95, but it is not consistent 
with the standard-theoretical prediction of Ay = 1. According to ([7]), we have 1/2 < p < 1 
in three dimensions. A reasonable numerical evidence for this relation has been obtained 
in |18j from the susceptibility data, where the MC estimate p = 0.555(17) has been 
reported for the 3D XY (0(2)) model. It corresponds to A^ = 1.929(21) according to ([7]). 



3 Critical point singularities in the n— vector models 

In vicinity of the phase transition point, various quantities have often power-law singular- 
ities, which are described by the critical exponents. For three-dimensional systems, exact 
results are difficult to obtain, and one usually relies on approximate methods. A review of 
standard perturbative RG methods applied here can be found, e. g., in [20" 1 |2~T | l22l |23~1 l24j . 
The results of the Borel-resummation of the perturbation series for the critical exponents, 
obtained within this approach, are reported in [25]. We will further focus on the Monte 
Carlo testing of the theoretical predictions for the 3D Ising model, which is a particular 
case of n = 1. The critical exponent rj, describing the ~ k~ 2+v singularity of the critical 
correlation function G*(k), is of particular interest here. According to [25] . the most ac- 
curate theoretical value r\ = 0.0335 ± 0.0025 is obtained from the series at fixed dimension 
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d = 3. The results of the resummation of the e-expansion in [25] are rj = 0.036 ±0.005 and 
T) = 0.0365 ± 0.0050. If all these estimates are correct within the error bars, then we have 
0.0315 < rj < 0.036. This estimation fairly well agrees with the value rj = 0.0366(8) ex- 
tracted from the finite-size scaling analysis of the MC data within L € [10, 128], reported 
in We will present the MC results for much larger lattices up to L = 1536 to test the 
agreement more precisely. The other relevant here exponents are the correlation length 
exponent v and correction-to-scaling exponent u. The widely accepted values for the 3D 
Ising model are v ~ 0.63 and ui ~ 0.8. They are in agreement, e. g., with the estimates 
v = 0.6304 ± 0.0013 and u = 0.799 ± 0.011 reported in [25]. 

As discussed in Sec. [2] the alternative theoretical approach, called the GFD theory, 
provides promising results concerning the Goldstone mode singularities. Therefore, it is 
interesting to verify the predictions of this theory also for the critical point singularities. 
In [16], a set of possible values of the exact critical exponents for the (p 4 model in two 
(d = 2) and three (d = 3) dimensions has been proposed: 

_ d + 2 j + Am _ 2(1 + m)+j 

7 ~ d(l + m + j) - 2j ' V ~ d(l + m + j) - 2j ' U 

where 7 is the susceptibility exponent, m may have a natural value starting with 1, and 
j is an integer equal or larger than —m. Other critical exponents can be calculated from 
these ones, using the known scaling relations. These values agree with the known exact 
solutions of the two-dimensional Ising model (m = 3, j = 0) and of the spherical model 
(J/m — > 00). A prediction has been made also for the 3D Ising model, i. e., 7 = 5/4 and 
v = 2/3, corresponding to m = 3 and j = 0, as in the two-dimensional case. This value of v 
is consistent with the logarithmic singularity of specific heat (according to a = 2 — du = 0) 
proposed earlier by Tseskis [27J. The exponents 7 = 5/4 and v = 2/3 have been later 
conjectured for the 3D Ising model by Zhang [28J. The critical exponents 7 = 9/8 and 
v = 5/8, calculated for the liquid-gas system by Bondarev [29], [30] , also correspond to ([5]). 
In this case m = 2 and j = —1 hold. As explained in [16] ; the equations ([8]) are meaningful 
for positive integer n, and we can have j = j{n) and m = m(n) in the case where the order 
parameter is an n-component vector. The spatial dimensionality d might be considered 
as a continuous parameter in ([8]) within 2<<i<4. At n = 1, the second-order phase 
transition at a finite temperature can be expected also for d < 2. However, according to 
the discussion in [31], an analytic extension of ([8]) to this region, probably, is only formal 
and does not correspond to any real system (fractal lattice). 

According to [16] . corrections to scaling can be represented by an expansion of cor- 
rection factor (amplitude) in integer powers of t 2v ~" f and t 2n< ~ du at t — > 0, where t is the 
reduced temperature. This expansion is simplified, since (27 — dv)/(2v — 7) is an integer 
number according to ([8]), and 2 — 7/1/ = ij holds according to the known scaling rela- 
tion. Hence, we obtain the expansion in powers of t 9 , where 9 = rju. In other words, the 
correction-to-scaling exponent in the expansions at t — > is 9 = r/i/, and the correspond- 
ing correction-to-scaling exponent in the finite-size-scaling analysis is uj = 0/v = rj, if 
the first expansion coefficient is nonzero. Allowing that some of the expansion coefficients 
are zero, we can have 9 = lr]v, where I is a positive integer. However, our actual MC 
data support the most natural choice of I = 1 for the 3D Ising model. It implies that the 
expansion coefficients do not vanish in general, but the nontrivial correction terms vanish 
only in some special cases, as the 2D Ising model or the spherical model, where u = 1 
and /or 9 = 1 . According to the numerical transfer matrix calculations in [32j , a nontrivial 
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correction to finite-size scaling with the exponent u = rj, probably, exists in the two-point 
correlation function even in the 2D Ising model, although its amplitude is very small here. 

The discussed here GFD critical exponents r/ = 1/8 and v = 2/3 for n = 1 seem to 
be quite incompatible with the MC data of the 3D Ising model, if we assume the usual 
correction-to-scaling exponent oj ~ 0.8. However, it turns out that our data for large 
enough sizes are very well consistent even with the set of exponents rj = u) = 1/8 and 
v = 2/3. The disagreement of these exponents with those of the perturbative RG method 
can also be understood based on a critical analysis [33]. 

4 Monte Carlo simulation results for the 3D Ising model 

We have simulated the 3D Ising model on a simple cubic lattice (at h = 0), using the 
iterative method introduced in [33] to adjust the coupling (3 to the pseudo-critical coupling 
/3 C , which corresponds to certain value 1.6 of the quantity U = (m 4 ) /(m 2 ) 2 , where m is the 
magnetization per spin. The pseudo-critical coupling f3 c tends to the true critical coupling 
j3 c at L — > oo. By this method, we estimate the moments of energy e and magnetization 
m per spin, i. e., (e k m l ) for a set of (3 values fluctuating around /3 C , and then calculate 
these mean values and related quantities at any given f3 near (3 C by using the Taylor series 
expansion for hi((— e) k m l ). We have evaluated /3 C (L) for various lattice sizes L and then 
estimated /3 C by fitting these /3 C (L) data. Besides, we have calculated the susceptibility 
X = L 3 (m 2 ) and the derivative dQ/d(3, where Q = 1/U, at (3 = f3 c (L), as well as at certain 
estimated critical couplings f3 c , in order to make a finite-size scaling analysis of the critical 
exponents. 

This method has been tested and discussed in detail in [34] , giving also some data 
for L < 1024. Here we have extended the simulation results up to L = 1536. A parallel 
(OpenMP) implementation of the Wolff single cluster algorithm has been used in our 
simulations, combining it with a sophisticated shuffling scheme for generation of a high- 
quality pseudo-random numbers, as described in [34J. Several tests have been made 
in [34] to verify the quality of the pseudo-random numbers. Here we have tested and used 
also the lagged Fibonacci pseudo-random number generator (PRNG) with multiplication 
operation and the lags r = 24, s = 55 (see [35] )• We have improved it with the standard 
shuffling scheme |35j, using the shuffling box of the size 10 000 in this case. Although the 
lagged Fibonacci generator with additive operation produces certain correlations in the 
sequence of pseudo-random numbers, which can be well detected by the directed-random- 
walk test [36], this defect is practically not observed in the actual multiplicative case. Like 
in [34], we have verified it by performing such test with 10 12 trajectories. Good results 
have been obtained in this test both for the original PRNG and for the one improved by 
the shuffling. It indicates that the multiplicative lagged Fibonacci generator works fine in 
applications with cluster algorithms, as it also has been mentioned in [37j . 

The simulations of the two largest lattices with L = 1280 and L = 1536 have been 
performed with the shuffling scheme of [35j and partly (24 usable iterations from 56 or 54 
ones at L = 1280 and L = 1536, respectively) also with the improved lagged Fibonacci 
generator. For complete confidence, we have verified that the simulation results of the 
two generators agree within the statistical error bars, whereas the final simulated values 
have been obtained by summarizing the data from both of them. Note that one iteration 
included 5 280 000 MC (Wolff cluster algorithm) steps at L = 1280 and 6 720000 MC 
steps at L = 1536, the MC measurements being performed after each 160 and 192 steps, 
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Table 1: The values of /3 C , as well as x/L 2 and 10 3 dQ/df3 at (3 = j3 c depending on L. 



T 

J_l 


a 

Pc 


XI L 


iu OL^/op 


1536 


0.2216546081(114) 


1.1900(22) 


96.01(73) 


1280 


0.2216546524(136) 


1.2009(25) 


72.04(53) 


1024 


0.221654625(22) 


1.2046(28) 


50.57(45) 


864 


0.221654635(25) 


1.2165(21) 


38.71(25) 


768 


0.221654672(27) 


1.2212(20) 


31.89(20) 


640 


0.221654615(31) 


1.2281(17) 


23.95(12) 


512 


0.221654662(45) 


1.2367(16) 


16.785(77) 


432 


0.221654637(58) 


1.2450(18) 


12.907(59) 


384 


0.221654567(65) 


1.2480(16) 


10.627(50) 


320 


0.221654716(75) 


1.2578(16) 


7.967(36) 


256 


0.22165460(11) 


1.2656(15) 


5.577(24) 


216 


0.22165460(13) 


1.2726(12) 


4.288(14) 


192 


0.22165425(16) 


1.2734(14) 


3.533(14) 


160 


0.22165414(18) 


1.2818(11) 


2.6495(87) 


128 


0.22165430(20) 


1.2913(10) 


1.8643(49) 


108 


0.22165376(26) 


1.2969(10) 


1.4170(36) 


96 


0.22165369(32) 


1.3012(10) 


1.1796(28) 


80 


0.22165278(32) 


1.30659(74) 


0.8822(18) 


64 


0.22165159(52) 


1.31466(78) 


0.6192(11) 


54 


0.22164968(56) 


1.31916(76) 


0.47334(81) 


48 


0.22164790(69) 


1.32164(66) 


0.39331(63) 


40 


0.22164383(80) 


1.32562(64) 


0.29424(40) 


32 


0.22163444(98) 


1.32835(59) 


0.20586(25) 


27 


0.2216212(11) 


1.32829(52) 


0.15703(16) 


24 


0.2216125(12) 


1.33027(47) 


0.13076(13) 


20 


0.2215821(17) 


1.32860(42) 


0.097717(87) 


16 


0.2215235(18) 


1.32510(34) 


0.068538(48) 



respectively. 

Our simulation results for the pseudo-critical coupling /3 C , as well as for the correspond- 
ing values ofx/L 2 and dQ / dfi depending on the lattice size L are given in Tab. [TJ In Tab. [21 
the latter two quantities are given, calculated at f3 = 0.221654604 and ft = 0.221654615, 
corresponding to two different estimates of the critical coupling. Besides, we have calcu- 
lated these quantities also at f3 = 0.2216545 for L < 128 and have verified that the values 
agree within the error bars with those reported for this (3 in [26J. 

5 Estimation of the critical coupling 

We have evaluated the critical coupling (3 C by fitting the j3 c {L) data to the finite-size 
scaling ansatz 

ft c (L) ~{3 C + L- l l» (a + ai L~" + a 2 L- 2ui ) , (9) 
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Table 2: The values of x/L 2 and 1(T 3 dQ/d(3 at (3 = 0.221654604 and p = 0.221654615 
depending on L. 





P = 0.221654604 


P = 0.221654615 


L 


y/L 2 
XI *-* 


10~ 3 dO/dB 


y/L 2 
XI u 


\Q~ 3 dO/dB 


1536 


l.looz(oo J 


CiK CIA ( ftO"\ 

yo.y4(by j 


1 1 CiQCii 

l.iyoU(oo J 


Oft 1 9/ ftO"\ 

yo.iz(by j 


1280 


l.lo4y(oU J 


71 ROf AK\ 
1 1.0/(40 J 


l.looO^oU J 


71 79^/1 K\ 
I L.I Z(40 J 


1024 


i.iyyo(o/ ) 


cn A o( Ar\\ 
011.40(41) J 


1 909^^7^ 
l.ZUZo(o / ) 


0U.0o(4U J 


864 


1911 r\('i r >\ 

l.ZllU(OZ J 


oq ftQ/'99 N \ 
OO.Oo(ZZ ) 


1 91 '?n('Q9 s \ 
l.ZloU(OZ J 


OO.Db(ZZ J 


768 


1 91 1 (V Q1 
l.ZllU(ol J 


01 77/ 1 s ! 
oL.I I (loj 


1 91 9fi('Q9 s \ 
!.ZlZb(oZ J 


01. / y(io j 


640 


1 99ftQ^'97^ 
l.ZZbO(Z f J 


9Q O/l /i 
Zo.y4(lU J 


1 99Q^^'97^ 

l.ZzoU(z 1 ) 


io. yo(iu j 


512 


1 0*301 

l.zozl(z7 ) 


1ft 7£ft / ftO\ 

lb. /ob(by J 


1 9QQ^/'97^ 

l.zooU(z7 ) 


1 ft 7ft9 ^ ftfA 

lb. r bz(by J 


432 


1 9/L9Ql'9zL N l 
i.z^izyi z^± i 




l.Z^lOUl Z^± I 


1 9 QD1 t^9) 

iz.yui (oz j 


384 


1.2499(24) 


10.634(44) 


1.2505(24) 


10.636(44) 


320 


1.2534(22) 


7.955(33) 


1.2538(22) 


7.956(33) 


256 


1.2656(23) 


5.577(22) 


1.2659(23) 


5.577(22) 


216 


1.2727(21) 


4.288(13) 


1.2729(21) 


4.288(13) 


192 


1.2795(21) 


3.540(12) 


1.2797(21) 


3.541(12) 


160 


1.2879(17) 


2.6550(78) 


1.2880(17) 


2.6551(78) 


128 


1.2941(13) 


1.8660(42) 


1.2942(13) 


1.8661(42) 


108 


1.3028(13) 


1.4198(32) 


1.3029(13) 


1.4198(32) 


96 


1.3066(14) 


1.1817(25) 


1.3066(14) 


1.1817(25) 


80 


1.3147(11) 


0.8845(16) 


1.3147(11) 


0.8845(16) 


64 


1.3241(12) 


0.62115(96) 


1.3241(12) 


0.62116(96) 


54 


1.33098(85) 


0.47518(69) 


1.33101(85) 


0.47518(69) 


48 


1.33504(93) 


0.39503(54) 


1.33506(93) 


0.39503(54) 


40 


1.34177(77) 


0.29577(34) 


1.34179(77) 


0.29577(34) 


32 


1.34958(71) 


0.20728(21) 


1.34960(71) 


0.20728(21) 


27 


1.35514(60) 


0.15838(14) 


1.35515(60) 


0.15838(14) 


24 


1.35836(53) 


0.13192(11) 


1.35836(53) 


0.13192(11) 


20 


1.36471(61) 


0.098821(75) 


1.36471(61) 


0.098821(75) 


16 


1.37072(45) 


0.069508(41) 


1.37072(45) 


0.069508(41) 
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0.221654 
0.221653 
0.221652 



0.0005 0.001 -1/0.63 

.Li 

Figure 1: The pseudo-critical coupling f3 c vs L _1// ° 63 . The dashed line shows the fit to ([!]), 
including only the leading correction term to f3 c — f3 c with coefficient a\ and exponents 
v = 0.63 and to = 0.8, whereas the solid line is the fit with both correction terms included. 

0.221654 
0.221653 
0.221652 



0.0005 0.001 0.0015 -1.5 

Figure 2: The pseudo-critical coupling j3 c vs L~ 3 / 2 . The solid curve shows the fit to ([9]), 
including only the leading correction term to f3 c — f3 c with coefficient a\ and exponents 
v = 2/3 and u = 1/8. 

neglecting higher order corrections. The fits with the usual (RG) exponents v = 0.63 
and oj = 0.8 within the range L 6 [64, 1536] are shown in Fig. [TJ The dashed line shows 
the fit with only the leading correction term (the term with coefficient a\) to (3 C — f3 c 
included, whereas the solid line - the fit with both correction terms. These fits give 
fi c = 0.2216546234(99) and /3 C = 0.221654615(13), respectively. The x 2 / d -°- f - of both these 
fits is 1.08. The data for L < 64 are omitted here, since inclusion of these relatively small 
sizes only increases the systematic errors without an essential reduction of the statistical 
errors. 

Surprisingly, the j3 c data within L € [64, 1536] can be even better fit with the exponents 
v = 2/3 and ui = 1/8 of the GFD theory. The fit, with only the leading correction to 
/3 C - P c included, is shown in Fig. El It yields j3 c = 0.221654585(15). In this case x 2 / d -°- f - 
has the value 0.97, which is even smaller than for the fits with v = 0.63 and u> = 0.8. 
Inclusion of the second correction term changes the result only slightly, i. e., it gives 
(3 C = 0.221654588(47) with x 2 /d.o.f. = 1.03, and the fit curve lies almost on the top of 
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eff 




0.04 



0.035 



0.03 



0.45 



0.5 



0.55 



L 



-1/8 



Figure 3: The effective exponent r/ e fj vs L 1 / 8 at (3 = 0.221654604 (solid circles) and 
P = 0.221654615 (x). The solid line is the fit of circles to CLE} with rj = u = 1/8. 

that one shown in Fig. [2j The statistical error, however, is strongly increased in this case. 

Alternatively, the critical coupling j3 c can be determined from the susceptibility data, 
requiring the consistency with the finite-size scaling ansatz 



at P = p c . Here the effective exponent rj e g{L) is the local slope of the — In (x/L 2 ) vs InL 
plot, evaluated from the linear fit within [L/2, 2L\. We observe that the plot of the effective 
exponent within L £ [64, 768] (evaluated from the data within L € [32, 1536]) is optimally 
described by fill j) with ?7 = u; = l/8at/3~ 0.221654604, as shown in Fig. [3l From this 
we conclude that, if the exponents rj = cu = 1/8 are correct, then p c = 0.221654604(18). 
This is the estimate of P at which the fit to fjll j) with cj = 1/8 and rj as a fit parameter 
yields rj = 1/8 within the statistical error bars. The solid-line fit at P = 0.221654604 is 
really good, and the effective exponents can be quite well fit with rj = oj = 1/8 also at 
P = 0.221654615 (the data shown by "x"), which is the best p c value provided by the 
already considered estimation with the RG exponents. 

The value p c = 0.221654604(18) is our best estimate of the critical coupling in the case 
where we use the exponents of the GFD theory, since two corrections to scaling are in- 
cluded, giving smaller statistical error bars than for the other estimate p c = 0.221654588(47) 
with two corrections included. Our estimated value P c = 0.221654604(18) is sufficiently 
reasonable and plausible, since it agrees well within the error bars with all other estima- 
tions considered here, as well as with the most accurate values of p c provided by other 
authors, i. e., p c = 0.22165455(5) [38] and p c = 0.22165457(3) [39]. 

Concerning the standard RG exponents, the method of fitting rj e g data is less conve- 
nient, since it is sensitive to the precise value of rj, which is not known accurately enough 
in this approach. 




(10) 



or with the corresponding ansatz for the effective exponent 



Tfaff(L) ~?? + ciL u + c 2 L 



(11) 
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Figure 4: The effective exponent 7] c g vs L~ 0,8 at /3 = /3 C (L). The solid curve is the fit 
to (fTTj) with fixed uj = 0.8. The dotted lines indicate the range [0.0315, 0.036] of rj values 
consistent with the set of perturbative RG estimations in |25| . 

6 Finite— size scaling analysis of the critical exponents 
6.1 Tests of the RG exponents 

In order to test the results of the perturbative RG theory, we have estimated the critical 
exponent r\ by fitting the susceptibility data to (fTUj) and also by fitting the rj e s data to (fTTj) 
at fixed RG exponent oj = 0.8. Note that (flOl) and (fTT|) are valid at f3 = f3 c , as well as at 
j3 = (3 C (L). The fits have been made at /3 = /3 C (L) and at the estimated in Sec. [5] values of 
the critical coupling, obtained by using the RG exponents. 

We start with the fitting of the effective exponent, since this method can be better 
controlled to see which fits are most appropriate. The plot of 77 e ff vs with uj = 0.8, 
evaluated at (3 = /3 C (L), is shown in Fig. H] We see that the data are well described by a 
quadratic curve for L > 108, therefore one has to include both correction terms in (jlip . 
The fit yields r] = 0.0380(26). This value is slightly larger than the upper limit 0.036 
of the interval 0.0315 < r/ < 0.036 for the values, which are consistent with the set of 
perturbative RG estimations of [25] discussed in Sec. El However, the discrepancy is 
within the error bars. Somewhat larger deviations above 0.036 are provided by the fits 
within L > 64, shown in Fig. GJ at approximately estimated values f3 c = 0.2216546234(99) 
and p c = 0.221654615(13) of the critical coupling. These fits yield 77 = 0.0381(18) and 
77 = 0.0397(28), respectively. The indicated here standard errors are calculated by the 
jackknife method, taking into account the statistical fluctuations in the estimated f3 c value. 
Namely, the standard error a is evaluated as a = JYUj A^, where Ay is the shift in the 
estimated r/ value, taking into account also the shift in the fitted (3 C value, when the j'th 
simulation block (iteration) for the ith lattice size is omitted. The statistical correlations 
are such that the total standard error is smaller than in the case where there are no 
correlations between the fluctuations in the estimated j3 c value and in the r] value obtained 
at a given j3. 

We have also evaluated 77 from the fits to §W\\ at the pseudo-critical coupling /3 C (L), as 
well as at the two estimates /3 C = 0.2216546234(99) and f3 c = 0.221654615(13) of the critical 
coupling. The results of fits within L G [L m i n , 1536], assuming u = 0.8, are collected 
in Tab. El The best estimates in these three cases are assumed to be 77 = 0.0380(23), 
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Figure 5: The effective exponent 7] cS vs L~°- 8 at f3 c = 0.2216546234(99) (diamonds) 
and p e = 0.221654615(13) (circles). The curves are the fits to (03) with fixed u = 0.8. 
The dotted lines indicate the range [0.0315, 0.036] of rj values consistent with the set of 
perturbative RG estimations in [25J. 



Table 3: The critical exponent ry evaluated from the fits to (jlOh within L E [£ mm , 1536] 
with fixed uj = 0.8. The estimates from the data at f3 = /3 C (L) and at /3 = /3 C with f3 c = 
0.2216546234(99) and p c = 0.221654615(13) are denoted by rfr, m and 773, respectively. 



^min 


m 


m 


m 


32 


0.0345(12) 


0.0363(14) 


0.0379(25) 


40 


0.0353(14) 


0.0364(15) 


0.0384(27) 


48 


0.0359(17) 


0.0369(17) 


0.0394(29) 


54 


0.0369(19) 


0.0376(20) 


0.0402(28) 


64 


0.0380(23) 


0.0376(30) 


0.0405(25) 


80 


0.0377(37) 


0.0366(38) 


0.0398(32) 
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Figure 6: The effective exponent r/ e fj vs L 1//8 at /3 = /3 C (L). The solid line is the fit with 
rj = uj = 1/8. 

77 = 0.0376(20) and 77 = 0.0405(25), since they perfectly agree with the corresponding 
estimates at a larger L m - m , but have smaller statistical errors. Note that the decrease of 
the standard error for 773 in Tab. [3] at L m \ a = 64 as compared to L m \ n = 48, 54 is the effect 
of the already mentioned statistical correlations. These best 77 values are very similar to 
the ones obtained before from the effective exponents. 

Remarkable is the fact that a better estimation of the critical coupling, including two 
corrections to scaling instead of only one correction, leads to a worse agreement with the 
results of the perturbative RG theory. In particular, the value 77 = 0.0405(25) deviates 
above 0.036 (the upper limit of the best RG estimate 77 = 0.0335 ± 0.0025 of [25]) by l.Sa. 
The observed here deviations suggest that 77, probably, is larger than normally expected 
from the perturbative RG theory. The statistical errors, however, do not allow to make a 
strict conclusion. 

6.2 Fits with the exponents of the GFD theory 

As already discussed in Sec. El the pseudo-critical coupling f3 c (L) can be even better fit 
with the exponents v = 2/3 and oj = 1/8 of the GFD theory than with those of the 
perturbative RG theory. We have also verified that the effective exponent r) e g(L) is very 
well described by the ansatz (fTTj) with 77 = uj = 1/8 (see Fig. [3]) at a certain /3 corresponding 
to a reasonable estimate of the critical coupling j3 c = 0.221654604(18). Also at the pseudo- 
critical coupling /3 = /3 C (L), the effective exponent r] e g(L) within L € [108,768] (extracted 
from the susceptibility data within L 6 [54, 1536]) can be well fit to (llip with fixed GFD 
exponents 77 = oj = 1/8, as shown in Fig. [6j If the exponent 77 is considered as a fit 
parameter at fixed oj = 1/8, then this fit yields 77 = 0.079(80) in agreement with the 
expected value 1/8, although the statistical error is rather large in this case. Sufficiently 
accurate data for even larger L values would be very helpful here to perform a more precise 
testing of the agreement. 

We have verified that not only the effective exponent 77^, but also the susceptibility 
data can be well described by the exponents 77 = oj = 1/8. In particular, x 2 /d.o.f of the 
fit with 77 = uj = 1/8 is 0.88 for the susceptibility data at (3 = (3 C (L) within L 6 [64, 1536]. 
At j3 = 0.221654604 ~ /3 C , the data can be well fit with x 2 /d.o.f = 1.01 over a wider range 
of sizes L G [48,1536]. 
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We have fitted also the 8Q/df3 data to the finite-size scaling ansatz 

|| ~ L x l v (a + A X L-" + A 2 L~ 2 ") . (12) 

These data can be well fit with the RG exponents, as well as with those of the GFD 
theory. In the latter case (y = 2/3, u = 1/8), the data within L € [32, 1536] are fit with 
X 2 /d.o.f = 0.88 at /3 = /3 C (L) and with x 2 /d-o.f = 0.93 at /3 = 0.221654604 ~ /3 C . 

In fact, reasonable fits with the exponents r\ = 1/8 and v = 2/3 are possible because 
of the value 1/8 of the correction-to-scaling exponent u. A remarkably larger value, such 
as u) ~ 0.8, would make such fits not good. According to (jlip . the plot of the effective 
exponent rj e g vs L~ u has to be almost linear at L — > oo. If to is as large as 0.8, then L~ 2uJ 
is already quite small for L > 108, so that a good linearity of the fit curve in Fig. [J] can 
be expected. However, the plot looks rather nonlinear. Moreover, it can be much better 
approximated by a straight line (within L € [108,768]) in the L _1//8 scale than in the 
^-o.s sca i e; as [i [ s dearly seen comparing the plots in Figs. 0]and[6l According to this, 
it seems, in fact, very plausible that u could be about 1/8. 



7 Conclusions 

The analysis of the MC data for the 0(2) and 0(4) models below the critical point sup- 
ports the recently published theoretical results [15], predicting that the Goldstone mode 
singularities in the 0{n) models are described by nontrivial exponents, as discussed in 
Sec. [2j Therefore, it has been important to verify the earlier predictions of this approach, 
called the GFD theory, concerning the critical point singularities. For this purpose, we 
have performed MC simulations of the 3D Ising model for very large lattices with linear 
size up to L = 1536, using a parallel implementation of the Wolff single cluster algorithm. 
The finite-size scaling analysis shows that the actual data can reasonably well interpreted 
with the usual critical exponents rj ~ 0.0335, v ~ 0.63 and ui ~ 0.8 of the perturbative RG 
theory, and can also be well described by a different set of critical exponents, 77 = ui = 1/8 
and v = 2/3, found within the GFD theory. The validity of the fits with the latter set 
of exponents depends on whether the correction-to-scaling exponent u; can, indeed, be as 
small as 1/8. We have demonstrated in Sec. 16.21 that this is a very plausible scenario. 

We have performed a high accuracy estimation of the critical coupling f3 c , combined 
with a sophisticated fitting of the data, including the leading, as well as sub-leading, 
corrections to scaling, to distinguish between the two sets of the critical exponents. We 
have found that a self-consistent estimation of f3 c and 77, assuming the RG exponents 
v ~ 0.63 and u = 0.8, gives somewhat larger values of r] (e. g., 77 = 0.0397(28) and 
77 = 0.0405(25) when two corrections to scaling are included) than the value 0.0335±0.0025 
expected from the most accurate resummation of the RG perturbative series [25] . These 
discrepancies are not large enough to make a strict conclusion about an inconsistency. On 
the other hand, the observed gradual increase in the effective exponent r] e f[(L) with growing 
system size L is just expected if 77 = a; = 1/8 and v = 2/3 are the correct asymptotic 
exponents. Therefore, it would be very interesting and important to see whether such an 
inrease to even larger t^s values is supported by the data for even larger lattice sizes. 
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